Exploring dynamic localization with a Bose-Einstein condensate 
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We report on the experimental observation of dynamic localization of a Bose-Einstein condensate 
in a shaken optical lattice, both for sinusoidal and square-wave forcing. The formulation of this 
effect in terms of a quasienergy band collapse, backed by the excellent agreement of the observed 
collapse points with the theoretical predictions, suggests the feasibility of systematic quasienergy 
band engineering. 
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I. INTRODUCTION 

The seminal work by Dunlap and Kenkre [l| has led 
to the recognition that an external uniform ac force can 
control the spreading of the wave packet of a particle 
moving in a spatially periodic potential. In particular, 
an initially localized wave packet of a monochromati- 
cally driven particle on a nearest-neighbor tight-binding 
lattice remains perpetually localized for certain values 
of the driving amplitude This phenomenon, termed 
"dynamic localization" , is a quantum mechanical mani- 
festation of the fact that a time-periodic force sometimes 
can stabilize a system, as shown in classical mechanics 
by the example of the driven, inverted pendulum 
Dynamic localization is closely related to the coherent 
destruction of tunneling experienced by a single parti- 
cle in a double-well potential under the influence of an 
ac force [1, S H @|, and to the modification of atomic 
(^-factors in oscillating magnetic fields 0, view 
of possible applications to electrons in terahertz-driven 
semiconductor superlattices, it has been shown theoret- 
ically that dynamic localization survives in the presence 
of Coulomb interactions [9]. However, experiments with 
such devices are difficult to perform, and have to cope 
with a host of competing effects While there ex- 

ist clean visualizations of dynamic localization in optical 
analogs of driven quantum systems, realized by means 
of curved waveguide arrays for light [ll], [l^] , the obser- 
vation of dynamic localization of matter waves has long 
remained a challenge. 

The situation changed with the availability of ultra- 
cold atoms in optical potentials. In 1998 Madison et 
at obtained evidence for band narrowing with cold 
sodium atoms in a phase-modulated optical lattice, in 
good agreement with theoretical calculations going be- 
yond both the tight-binding and the single-band approx- 
imation , 13J . More recently, the basic mechanism respon- 
sible for dynamic localization, an effective rescaling of the 
hopping matrix element induced by the ac force, has been 



observed for single-particle tunneling in strongly driven 
double- well potentials fT3| . It had also been pointed out 
that this mechanism remains effective even for an inter- 
acting Bose-Einstein condensate, at least for sufficiently 
high driving frequencies fl^ . This has led to the experi- 
mental demonstration of dynamic control of matter- wave 
tunneling in an optical lattice [l^, and to the observa- 
tion of photon-assisted tunneling of a condensate 
Since the reduction of interwell tunneling increases the 
relative importance of the particles' repulsion, even the 
superfluid-to-Mott insulator transition of ultracold atoms 
in an optical lattice can be induced by shaking the lat- 
tice in a time-periodic manner, as has now been shown 
by Zenesini et al. [l^ 

In the present paper we take up this line of investiga- 
tion and explore dynamic localization of a dilute Bose- 
Einstein condensate in a time-periodically shifted opti- 
cal lattice in more detail. We proceed as follows: In 
Sec. [H] we cast the concept of single-particle dynamic lo- 
calization into a form that lends itself in a particularly 
transparent manner to further gen eralizations , relying on 
the idea of quasienergy bands . We then apply these 
general considerations in Sec. IIIII to optical cosine lat- 
tices, taking into account all couplings that are omitted 
in the nearest-neighbor approximation. In Sec. llVl we re- 
port our experimental results, achieved with condensates 
of ^^Rb. Besides sinusoidal forcing, we also consider the 
case of square-wave forcing, which is known to produce 
exact single-particle dynamic localization for any form of 
the energy dispersion [l^, [2l|. The final Sec. |V] contains 
our conclusions. 



II. THE PRINCIPLE UNDERLYING DYNAMIC 
LOCALIZATION 

We consider a time-dependent Hamiltonian 

H{t)^Ho+Hi{t) , (1) 
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where 



(2) 



describes a single particle with mass m moving in a one- 
dimensional periodic potential V{x) = V{x + d) with 
lattice constant d, and 



Hi{t) = ~F{t)x 



(3) 



introduces an external force F(t), mediated by the po- 
sition operator a; in a manner analogous to electromag- 
netic forces on atoms within the dipole approximation. 
The energy eigenf unctions 'fn,k{x) of i?o are Bloch states 
with band index n and wave number fc, 



(4) 



with functions Vn_k{x) — Vn,k{x + d) sharing the spatial 
lattice period d. The solutions to the eigenvalue equation 



Hoipn,k{x) = En{k)ipri,k{x) 



(5) 



provide the dispersion relations En(k) for the various en- 
ergy bands. It is of interest to observe that the external 
forcing does not truly break translational symmetry: In- 
troducing the unitary transfomation 



U 



cxp^— — ^ drF(T)a;^ 



the transformed wave functions 'ip{x,t) 
governed by the Hamiltonian 



(6) 

U'ilj{x,t) are 



P 



dTF{T)] +V{x), (7) 



which still remains invariant under spatial translations 
by d. ^ 

In the following, we assume band gaps so wide that, 
despite the external force, the dynamics can be restricted 
to the lowest band. Dropping the band index, we expand 
the Bloch waves of that band with respect to the corre- 
sponding Wannier states [13, HI] , 



fkix) 



wi{x)e'' 



(8) 



where wi{x) = wo(x — £d) denotes a Wannier function 
centered around the i-th lattice site. We take the poten- 
tial to be symmetric, V{x) = V{—x), so that wqIx) can 
be chosen real and symmetric 23i]. Defining the time- 
dependent wave number 



1 /•* 

qk{t)=k+- / dTF(r) 



(9) 



obeying hq^.{t) — F{t) and containing a still unspecified 
lower integration bound the wave functions 

=^«;,(x)e'«'=(*)"exp('-i fdTE{qk{T)) 

(10) 



then are "weak" solutions to the time-dependent 
Schrodinger equation with the full Hamiltonian ([T]), in 
the following sense: By construction, one has 



d 

ih—'tPkix,t) ^ E{qkit))tPkit) 



(11) 



-F{t)Y, w, ix)e"^' exp i ^*dT (g,. (r)) 



If we stipulate that the lattice site labeled £ — is situ- 
ated at a; = 0, implying 



(0|x|0} = dxx\wo{x)\^ = 



then one has 



{£\x\m) 



£d 




if £ : 

else 



TO 



(12) 



(13) 



where \£) is the Dirac ket corresponding to wg{x). Note 
that this is an identity, which holds exactly even for shal- 
low lattices. Hence, it follows that 

{£\ih^^Mx,t)) = {£\E{qk{t)) -F{t)x\Mx,t)) 

= {£\Ho + Hiit)\4,k{x,t)) (14) 

for each £. 

The above wave functions ([TU)) can be regarded as 
Houston states [24|, also known as "accelerated Bloch 
states", and apply to any type of uniform forcing F{t), 
provided the single-band approximation remains viable. 
If we now require that the forcing be periodic with pe- 
riod T and zero average, so that F{t) — F{t -f T) and 
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(15) 



a further step can be made: In this case the transformed 
Hamiltonian ([7]) is periodic in both space and time, thus 
giving rise to spatio-temporal Bloch waves 

■0fe(a;, t) — Vk{x, t) exp {ikx — ie{k)t/fij (16) 

with quasimomenta hk, quasienergies e(fc), and time- 
dependent Bloch functions obeying v{x, t) — v{x + d, t) = 
v{x,t + T). 

When deriving these states within the single-band ap- 
proximation from the Houston states (|10p . a slight sub- 
tlety comes into play: The lower integration bound to in 
the definition ^ of qk (t) effectuates a shift of the wave 
number k. In order to avoid this shift, and thus to make 
sure that the wave number k labeling a spatio-temporal 
Bloch wave is the same as the one which labels the Bloch 
state continuously connected to it when the driving force 
vanishes, this lower bound to has to be chosen such that 
the integrated force also has zero average, requiring 



1 fdTF{T)^o. 

^ Jo Jto 



(17) 
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With this specification of qk{t), we set 



£{k)^^J^ dtE{qk{t)) , 
and write the Houston states (IIOD in the form 



i'k{x,t) 



Mfc(a;,i)exp( --e{k)t 



(18) 



(19) 



This formal-looking manipulation is of key significance: 
After the average phase growth has been factored out, the 
remaining functions Uk{x, t) acquire the temporal period- 
icity of the force, 



Uk{x,t) 
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TABLE I: Hopping matrix elements ce = {0\Ho\£)/Er for the 
lowest band of an optical cosine lattice (|25|l . as functions of 
the scaled lattice depth Vo/Ei. Observe that the sign of ce 
alternates with I. 



X exp 



1 

^ JO 

Uk{x,t + T) 



At 



{E{qu{T)) ~ e{k) 



(20) 



Thus, these wave functions (|19p are Floquet states [3, 
[l9l] , from which the spatio-temporal Bloch waves ((T6)l are 
obtained by applying the unitary transformation ([6|). At 
times t = io + sT (with integer s) they coincide with the 
Bloch waves ([8]), apart from a phase factor. Any single- 
band wave packet ip{x,t) can be expanded with respect 
to these basis states, 

1/2 „7r/d 



-TT/d 



dk a{k)uk{x,t) exp[ ——s{k)t 



(21) 

with some density a(k) which is time-independent de- 
spite the presence of the ac force. The T-periodic nature 
of the functions Uk{x,t) directly reflects the short-term 
response of such a packet to the force F{t), whereas the 
quasienergy dispersion e(fc) governs its long-time aver- 
aged evolution in the same manner as the original energy 
dispersion E{k) governs its evolution when F(t) = 0. For 
example, the average velocity Vkg of a wave packet ini- 
tially centered around some wave number fcg is given, in 
the presence of the forcing, by 



1 dejk) 
h dk 



(22) 



k — kn 



Moreover, the representation ((2T|) directly leads to a quite 
general and intuitive condition for dynamic localization: 
Exact dynamic localization obviously occurs when the 
quasienergy band collapses, so that all quasienergies co- 
incide, e{k) — Eq for all k, implying that the effective 
mass of a driven Bloch particle becomes infinite. All 
Floquet states comprising the wave packet (PT|) then ac- 
quire exactly the same phase in the course of one driving 
period, since 



ip{x, t) — exp 



1/2 „7r/d 



7r/d 



dfc a{k)uk{x, t) , 
(23) 



so that the packet can neither move (apart from its resid- 
ual T-periodic motion) nor spread, but simply repro- 
duces itself T-periodically, regardless of the form of the 
density a{k). When this condition of zero width of the 
quasienergy band cannot be met exactly, but there still is 
significant band narrowing, wave packet spreading can at 
least be reduced substantially by the external force (19l] . 



III. APPLICATION TO OPTICAL LATTICES 

For a lattice with inversion symmetry, the original en- 
ergy dispersion takes the form 

oo 

E{k) ^ Eo + ^2{0\HQ\e)cos{lkd) , (24) 
1=1 

with coefficients given by matrix elements of iJo between 
Wannier states located (. sites apart from each other. In 
particular, for an optical cosine lattice of the form ^] 



V{x) — — cos{2ki^x) 



(25) 



where is the wave number of the lattice-generating 
laser radiation and the depth Vq is proportional to the 
intensity of that radiation [2^, the Wannier states can 
easily be computed numerically [2^ |23| ■ Table [J lists 
the resulting coupling coefficients (0|iJo|-^) for the lowest 
band as functions of the scaled lattice depth Vq/Ei-, with 



2m 



(26) 



denoting the single-photon recoil energy of an atom with 
mass m. The relative error committed when adopting 
the commonly used nearest-neighbor approximation, i.e., 
when neglecting all {0\Hq\£) with ^ > 2, is on the order 
of 10% for an optical lattice with Vo/Ey = 3, and reduces 
to 1% only when Vo/E, « 10 [27*1. 
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A. Sinusoidal forcing 

We now consider purely sinusoidal forcing, as specified 
by F{t) — Fq cos(ujt + ip) with some arbitrary phase ip. 
In this case, Eq. ^ together with the requirement ITTl) 
yields 



Fq 

gfc(t) = fc + — sin(tjt + cp) , 



implying 



- dt cos {qk{t)£d) = Joi^Ko) cos{£kd) 



(27) 



(28) 



where Jo(z) is a Bessel function of order zero, and we 
have introduced the scaled driving amplitude 



Ko = 



Fod 
huj 



(29) 



Hence, the quasienergy dispersion resulting from the 
Bloch band ([M)) under sinusoidal forcing reads 



e{k) = + X] 2(0 I I i) -hieKo) cos(^fcd) . (30) 

For the particular example of quasienergy bands originat- 
ing from the lowest Bloch band of an optical lattice (pS)) . 
characterized by the matrix elements collected in TablelH 
the quasienergy band widths are depicted in Fig. [1] as 
functions of Kq for some typical depths Vq. Due to the 
dominance of the nearest-neighbor hopping matrix ele- 
ment (0|^fo|l) there is strong band narrowing for values 
of Kq close to the zeros of the Bessel function Jq, but the 
nonvanishing longer-range hopping elements prevent the 
band from collapsing completely, as emphasized by the 
inset. Thus, one expects appreciable, but incomplete dy- 
namic localization in sinusoidally driven shallow optical 
lattices. 



B. Square- wave forcing 

For square-wave forcing of the form 

m = 



Fo , < t < r/2 

Fo , T/2<t<T, 



Eqs. dH) and HH) yield 

k + Foit-T/4)/h , 0<t<T/2 



k + Fo(3T/4 - t)/h , T/2<t<T 



(31) 



(32) 



to be continued T-pcriodically to all t. This gives 



- / dt cos {qk{t)ed) 



smc 



cosiikd) , (33) 
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FIG. 1: Width Asin of the lowest quasienergy band of an 
optical cosine lattice under sinusoidal driving as function of 
the dimensionless amplitude ((29}, for Vo/Ei = 2 (dots), 3 
(long dashes), 5 (short dashes), and 10 (full line). The inset 
quantifies the extent of band narrowing for values of .ft'o close 
to the first zero jo,i = 2.405 of the Bessel function Jq. 




with being defined according to Eq. ([29]), having set 
Lo = 2t:/T. Moreover, we write sinc(z) = sin(2;)/z. Cor- 
respondingly, the quasienergy dispersion becomes 



£(/c) ^ Eo + J22{0\Ho\£) sine 



cos{ekd) . (34) 



Since all sinc-functions adopt their zeros simultaneously, 
there is a total collapse of this quasienergy band when 
Kq = 2v with integer v = 1,2,3,... . Thus, we recover 
the known fact that there is exact dynamic localization, 
regardless of both the values of the hopping matrix el- 
ements and the form of the wave packet, for these par- 
ticular driving amplitudes 20, 21]. Figure [2] shows the 
widths of quasienergy bands originating from the lowest 
Bloch band of an optical lattice under square-wave forc- 
ing; the total band collapse dX Kq = 2 is clearly visible 
in the inset. 



IV. EXPERIMENTAL RESULTS 

Our experimental setup is described in detail in 
Ref. [iBl- Briefly, we adiabatically load Bose-Einstcin 
condensates consisting of about 6 x 10* atoms of *^Rb 
into the lowest band of a one-dimensional optical lat- 
tice. The lattice is generated by focusing two counter- 
propagating linearly polarized laser beams of wavelength 
A = 27r/fcL = 842 nm onto the condensate, resulting 
in the periodic potential ([25]) along the beam direction. 
Each beam passes through an acousto-optic modulator, 
allowing us to introduce a frequency difference lS.v(t) be- 
tween the beams that can be used to accelerate or shake 
the lattice. In the laboratory frame of reference, the con- 
densate then experiences a time-dependent potential 



Viab(a;,i) = y cos( 2fcL 



[ dTA//(T) 

Jtl 



(35) 



5 



LU 



1.8 2.0 2.2 




0.04 



0.02 



0.00 



FIG. 2: Width Asqw of the lowest quasienergy band of an 
optical cosine lattice under square-wave driving as function 
of the dimensionless amplitude (|29p . for Vo/-Er ~ 2 (dots), 3 
(long dashes), 5 (short dashes), and 10 (full line). The inset 
illustrates that an exact band collapse occurs when Ko is a 
nonzero integer multiple of 2. 



By means of a unitary transformation to the comov- 
ing frame of reference, the single-particle Hamiltonian 
is brought into the form 

H = f + ^cosi2k,.)+m^^.. (36) 
2m 2 2 at 

Hence, prescribing an oscillating frequency difference 

Aiy{t) = Ai^max sin(tjt) (37) 

with amplitude Az^max, one obtains sinusoidal forcing 
with strength Fq = mujXAi>yniix/'2, amounting to the 
scaled driving amplitude 



Ko 



2 

On the other hand, the triangular protocol 



2Ai/„ 
2Ai/„ 



<t/T 

.(l-t/T) 



< < < T/2 
T/2<t <T 



(38) 



(39) 



results in square-wave forcing with strength Fq — 
mAAfmax/T, giving 

Ko^l^. (40) 

We first study the in situ expansion of a condensate 
under both kinds of forcing [l6|. To this end, initially 
an anisotropic, elongated harmonic trap with trapping 
frequencies of 20 Hz (longitudinally) and 80 Hz (radially) 
was superimposed on the lattice. After switching off the 
laser beam effectuating the longitudinal confinement, the 
condensate was free to expand in the direction of the 
shaken lattice. After some time, it was illuminated by 
a resonant flash, the shadow cast by which was imaged 
onto a CCD chip. 



We work with lattices possessing a fairly large depth 
between Vq — QE^ and Vb — QE^, so that the effect 
of longer-range hopping is not discernible at the scale 
of our experimental accuracy, and we may employ the 
nearest-neighbor approximation, keeping only the hop- 
ping matrix element J — — (0|7Jo|l) connecting adjacent 
lattice sites. Thus, the unperturbed energy band is well 
described by the cosine dispersion relation 



E(k) = -2Jcos(fcd) , 



(41) 



assuming that interparticle interactions remain negligi- 
ble. Accordingly, the quasienergy band of the driven sys- 
tem is approximated by 



e{k) = — 2,/off cos(/c(i) 



(42) 



where the effective hopping matrix element Jeff is gven 
by 



Joff = J Jo(-f^o) 



(43) 



for sinusoidal forcing, as follows from Eq. pO]). whereas 
Eq. dsn gives 



Jeff — Jsinc(7ri^o/2) 



(44) 



for a square- wave drive. The measured expansion rate 
d(Jx/dt of the condensate width along the lattice direc- 
tion is then to a good approximation proportional to the 
hopping element Jeff which effectively describes nearest- 
neighbor tunneling when the forcing is present [T6| . This 
is a working example of our general philosophy: The 
time-averaged evolution of the driven system proceeds 
in close analogy to that of an undriven one, with e{k) 
replacing E(k). In practice, rather than extracting Jeff 
from the expansion data, for each experiment we also 
measure the bare expansion rate (which is proportional 
to J) in the undriven system, from which the normalized 
effective hopping element Jeff /J is then calculated. 

In Fig. Owe plot our results for both types of forcing 
versus the scaled driving amplitude A'o, obtained with 
a sinusoidal drive of frequency uj/2tt = 1.0 kHz, and a 
square-wave one with uj /2t: — 1.5 kHz. Evidently, the 
single-particle expectations are matched quite well, with 
the data obtained for square-wave driving scattering a 
bit more strongly around the prediction (|44p than those 
for sinusoidal driving around the graph of Eq. (j43p . 

An alternative method for determining the points of 
maximum band collapse is to measure the phase coher- 
ence of the condensate in the shaken lattice . Switch- 
ing off the confining potential and letting the conden- 
sate fall under gravity for 20 ms, we obtain an interfer- 
ence pattern whose visibility reflects the condensate co- 
herence. Recording this visibility as a function of time, 
we extract the decay time constant Tdeph of the result- 
ing near-exponential function. In general, we find de- 
phasing times on the order of 100 ms in the presence of 
even strong sinusoidal driving, while they are somewhat 
shorter for the square-wave drive, as shown in Fig. 2) In 



FIG. 3: (Color online) Absolute value of the effective hop- 
ping matrix element Jctt for a condensate in an optical lat- 
tice with depth Vo/Ey = 6 under the influence of sinusoidal 
(open boxes) and square-wave (triangles) driving, normalized 
to the "bare" hopping element J which governs the conden- 
sate spreading in the undriven lattice. The lines correspond 
to the approximations (|43p and (|44|) . The driving frequencies 
uj/2n were 1.0 kHz for sinusoidal driving, and 1.5 kHz for the 
square- wave case. 

the immediate vicinity of the band collapse points, how- 
ever, the dephasing times are strongly reduced: When 
the quasienergy band width approaches zero, the indi- 
vidual lattice sites are effectively decoupled, so that the 
local phases evolve independently due to interatomic col- 
lisions, resulting in a rapid dephasing of the array. This 
effect leads to sharply pronounced dips in plots of Tdcph 
vs. Kq, allowing us to confirm the theoretically predicted 
first collapse points to fairly good accuracy: Kq = 2.4 for 
sinusoidal forcing, whereas Kq = 2.0 for the square-wave 
drive. 

While in the present work we are mainly concerned 
with identifying the conditions for partial or complete 
dynamic localization, useful insight can also be obtained 
by studying the time evolution of the width of the driven 
matter-wave packet away from the collapse points. In 
general, if the extension ax of a Floquet wave packet (|2T|) 
narrowly centered in fc-space around some ko is given by 
ax = a at time t = 0, that wave packet spreads in time 
according to 

~ 7a(fco)i, (45) 
with asymptotic, large-i expansion rate 

,M . ^ (46) 
na 

determined by the second derivative of the quasienergy 
dispersion at ko. According to Eqs. ([50]) and ([M]) . here 
the hopping elements {0\Ho\i) enter with weights pro- 
portional to each depending in a different manner on 
the driving amplitude. Hence, precise measurements of 
the expansion rate could yield information both on the 
state of the driven packet, i.e., on the wave number fco, 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Ko 

FIG. 4: (Color online) Dephasing times for a condensate in 
an optical lattice with depth Vo/Er —9 under the influence of 
square-wave (triangles; left scale) and sinusoidal (boxes; right 
scale) driving. 

and on the ncxt-to-nearest neighbor couplings. Further- 
more, while we have worked with condensate densities 
so low that interparticle interactions can be ignored and 
the single-particle picture remains applicable, it would be 
interesting to identify signatures of interparticle interac- 
tions in dynamic localization. 

V. CONCLUSIONS AND OUTLOOK 

In this work we have pointed out that wave packet 
motion and spreading in time-periodically forced spa- 
tially periodic potentials are governed by quasienergy 
bands e{k) in a manner which exactly parallels the de- 
scription of quantum dynamics in unforced lattices in 
terms of their energy bands E{k), and have experimen- 
tally demonstrated dynamic localization, i.e., "freezing" 
of the dynamics due to a quasienergy band collapse, with 
dilute Bose-Einstein condensates in optical lattices sub- 
jected to either sinusoidal or square- wave driving. Seen 
from a conceptual viewpoint, this experimental confir- 
mation of the ideas developed in the context of dynamic 
localization can be regarded as a proof of principle for the 
notion of quasienergy band engineering: Different types 
of forcing, applied to the same energy band, can lead to 
substantially different quasienergy dispersion relations, 
and hence can be employed for realizing band structures 
which even might not have an analog in traditional solid- 
state systems. 

In order to make contact with previous theoretical 
works [l|, [l^, HO, [m , we have restricted our investiga- 
tion here to forces with parameters which do not lead to 
significant coupling of several energy bands. In this case 
the construction of the quasienergy dispersion merely in- 
volves taking the average (HH]) and thus is an easy ex- 
ercise, since one is essentially dealing with just a free 
Bloch particle. However, such quasienergy bands, reflect- 
ing the existence of spatio-temporal Bloch waves (|16l) . 
originate from nothing more than the simultaneous pres- 
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ence of a spatially periodic potential and a temporally 
periodic force. Accordingly, they also exist when sev- 
eral unperturbed energy bands are coupled by the force; 
the quasienergy bands then exhibit fairly nontrivial fea- 
tures [8;]- The exploration of these further options for 
band engineering which result from deliberate interband 
coupling, together with the inclusion of interparticle in- 
teraction, should open up further interesting avenues. 
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